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Executive Summary of the Research 

Research extended the SRG systems design approach to higher levels of structural hierarchy, while 
incorporating dynamic character in the form of adaptive system concepts via a synergistic integration of 
materials science and applied mechanics. The research addressed both the methodology of rational design at 
new levels of materials complexity, and the fundamental principles of specific materials behaviors necessary to 
support quantitative design. The research addressed adaptive multilevel microstructures in two broad areas for 
applications of advanced alloys: 

(a) Gradient and Layered Systems with application in high performance gear technology and advanced dual¬ 
hardness armor, achieving novel strength/toughness combinations through control of coherent carbide 
precipitation in processing and tuned-stability crack tip transformation plasticity in service. 

(b) Biomimetic Laminate Systems combining an oxidation resistant and high-temperature creep resistant matrix 
alloy with a thermodynamically compatible ductile reinforcement consisting of a shape memory alloy imparting 
biologically-inspired ductile-crack-bridge toughening via thermoelastic martensite pseudoelasticity at low 
temperatures, and shape-memory-induced crack clamping for diffusional damage healing at high temperatures. 

The primary objectives of the research effort were the methodology and fundamental principles for rational 
design of complex materials. This included the principles of microstructural evolution during coherent carbide 
precipitation for efficient strengthening without embrittlement, the control of metastable austenite precipitation 
for optimal-stability transformation toughening, the behavior of pseudoelastic crack-bridge toughening, and the 
shape-memory-based healing of internal composite damage. Research emphasized the dynamics of 
microstructural evolution during processing and in service, and the design methodology for control of 
multilevel dynamic hierarchical microstructures as programmed adaptive systems. 

Specific property targets that guided the conceptual design of prototype materials for evaluation were: 1) A 
secondary hardening gear steel with a case hardness of Rc70 while maintaining the core of Kic > 60 MPa m at 
Rc 50 hardness up to 400 °C. 2) A dual hardness armor system with a front face hardness of Rc 65-70, and a 
transformation toughened back plate of Rc 55-57 and Kic > 80 MPa m’^^. 3) A shape-memory reinforced 
superalloy laminate system with a 10^ hr creep strength of 100 MPa at an operating temperature 600-700 °C. 

Exploratory research also addressed the conceptual design of multifunctional metal/ceramic hybrid armor 
systems, and transformation-plasticity based matrix alloys for shearable tungsten heavy alloy composites. 
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Background 

Materials as Hierarchical Systems 

The late C. S. Smith [1] expounded a "systems" view of materials as hierarchical structures with properties 
governed by dynamic evolution in processing and in service. The Smith philosophy has been developed into a 
methodology for rational design of complex materials, pioneered by the multi-institutional Steel Research 
Group (SRG) program centered at Northwestern University [2-4], demonstrating successful application in high 
toughness ultrahigh-strength alloy steels [4,5]. This view of a secondary-hardening martensitic alloy steel as a 
system is summarized in Figure 1 and 2 [2], depicting the hierarchy of microstructural subsystems controlling 
the required set of properties for desired performance, and the sequential processing subsystems which in turn 
control structure. The system chart for high power density gear systems, Figure 2, considers not only the 
case/core gradient system of the gear steel itself, but the integrated design of compatible superhard ceramic 
coatings and their interfaces, ultimately incorporating redesign of gear shape to optimally exploit new material 
properties. Also represented are alternative processing routes to provide affordable manufacturing strategies for 
specific applications. Organized within this framework, basic research integrating materials science, applied 
mechanics and quantum physics has focused on key process/structure and structure/property relations, defining 
quantitative microstructural objectives expressible in thermodynamic terms (including scaling factors governing 
dynamic evolution), and the THERMOCALC thermochemical database and software system [6] has been 
applied as an integrative design tool predicting multi-component alloy composition capable of achieving desired 
multiphase microstructures under prescribed processing conditions. The successful application of this approach 
in ultrahigh strength steels has made possible an undergraduate course in Materials Design [3], and the 
extension of thermodynamics-based systems design to ceramics and polymers is now being undertaken within 

the Department of Materials Science and Engineering at Northwestern [7]. 

C. S. Smith also emphasized the strong parallels between engineering materials and the hierarchical structures 
of the biological world [1,8]. The discovery that relatively primitive biological systems such as viruses and 
bacteria exploit such sophisticated materials technology as martensitic transformation plasticity and the shape 
memory effect [9] reinforces this parallelism. In response, an active area of current investigation is the field of 
biomimetics [10,11] in which complex biological systems are investigated to incorporate biological principles 
in the rational design of more sophisticated engineering materials. This has inspired the notion of extending 
materials structure to both higher levels of hierarchy and more dynamic character to achieve adaptive or 
"smart" microstructures programmed to evolve along predetermined pathways in response to a given stimulus to 
achieve novel properties. Such programming "brings to life" Smith's seminal vision of dynamic hierarchical 
materials. 

The Steel Research Group 

The science-based approach to the computational systems design of materials is based on the pioneering work 
of the multi-institutional interdisciplinary Steel Research Group (SRG) program directed by G. B. Olson. 
Initiated in 1985 under NSF support as a collaboration of MIT, Harvard, and Brown with government 
laboratories and industry, the SRG reorganized in 1988 with Northwestern as lead institution, funded by 
mission-oriented agencies including ARO, ONR, AFOSR, NASA and DOE with approximately 20% industry 
support, and continuing NSF sUpercomputing support. Progress in the first 5 years is summarized in the 
Proceedings of the 34th Sagamore Army Materials Research Conference [12], organized around the SRG effort. 
The SRG design approach integrating processing/ structure/property/performance relations within a multilevel 
hierarchical system structure, and the computational design tools stemming from the research integrating 
quantum physics, applied mechanics and materials science are described in a recent invited overview [13] and 
millennial essay [13a] for Science magazine, attached as Appendix A. Broadening the approach to design of all 
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materials was the subject of a recent Gordon Research Conference [14] and a recent international symposium on 
materials design [14a]. 
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Figure 1. System structure of secondary hardening martensitic alloy steel. 
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Figure 2. System structure of high power-density gears. 
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The design methodology 

As discussed further in Appendix A, Figures 1 and 2 concisely summarize the systems approach unifying our 
design methodology. The chart denotes the hierarchy of microstructural subsystems underlying the set of 
material properties necessary for desired performance, and the sequential stages of processing which govern 
their dynamic evolution. This systems view identifies and prioritizes the key structure/property and 
process/structure relations for which computational models are needed to support predictive design. The 
resulting hierarchy of computational tools is summarized in Appendix A (insert, p. 997). Dynamic models 
predict evolution of microstructure in both processing and service. Working down from the highest structural 
level, models support (a) transformation design controlling the evolution of martensitic microstructures and 
their dynamic competition with alternative diffusional transformations, (b) micromechanics design via 
simulation of unit fracture processes in service employing kinetics-based constitutive laws (c) nanodesign 
integrating continuum mechanics of interfacial coherency and capillarity with multicomponent diffusional 
precipitation kinetics for efficient strength control, and (d) quantum design integrating frontier computational 
quantum mechanics in control of interfacial adhesion phenomena underlying higher level behaviors. Output of 
these tools defines thermodynamic design parameters governing critical time and length scales, allowing use of 
the THERMOCALC thermochemical software and database system as the principal design integration tool to 
satisfy multiscale design requirements. The efficacy of the approach has now been demonstrated in numerous 
SRG designs described in Appendix A and several review articles [15-18]. The Materials-by-Design doctrine 
has gained international recognition [19]. 

Achievements 

Education Initiatives 

As MSE Associate Chairman, Prof Olson led a committee which restructured the department curriculum, 
preparing students to take Materials Design in the Junior year. A Senior level Process Design course developed 
by D. L. Johnson featuring process optimization via statistical design of experiments, is followed by a 2-quarter 
experimental Senior Project. This structure allows students to perform theory-based conceptual designs in the 
Junior year and evaluate their prototype material as a Senior Project using experimental methods of the Process 
Design course. The new design-oriented curriculum has been credited for Northwestern’s program tying MIT in 
the most recent US News and World Report rankings. As a member of the Engineering First Committee 
restructuring the early undergraduate engineering curriculum, Olson brought computational materials design 
projects to a new Freshman-level 2 quarter Engineering Design and Communications course, building on 
experience with high school student summer projects. A freshman team successfully designed and implemented 
a Sn-based shape-memory-reinforced self healing biomimetic alloy composite. The growing design culture has 
brought student entrepreneurship to the undergraduate level. A team of freshman design students who 
performed a software project started a website design company during their freshman year, and the 
organizational meeting of an Enterprise Club drew over 100 interested engineering students. In the past year, 
the club organized an Enterprise course that drew 90 students, culminating in presentation of their business 
plans to investors. 

The SRG brought a close synergy of research and teaching with the development of Olson's MSc C90 Materials 
Design course at Northwestern [3]. Student team projects in which undergraduates employ THERMOCALC 
computational thermodynamics to integrate results of SRG graduate research have won national student 
materials design competitions [20-23] and have been featured on PBS-TV [24], several magazine articles [25- 
27], and a recent popular book on materials [28]. Projects have involved collaboration with design classes in 
other departments in the schools of engineering and management, and establishment of a similar materials 
design course in Sweden fostered through doctoral student exchange has allowed international undergraduate 




collaboration via internet [29]. Management school collaboration has facilitated student entrepreneurship 
through market analysis and business plans. After graduation, Dr. Charles Kuehmann founded QuesTek 
Innovations LLC offering eomputational alloy design services based on SRG technology, and new commercial 
alloys are undergoing field tests. In December 1999, QuesTek received the Technology of the Year Award from 
Industry Week magazine, in recognition of the new "materials by design" technology. 

Fundamental principles for design of complex materials 

1. Strengthening 

The first dynamic problem of interest in the systems of Figures 1 and 2 is the evolution of coherent carbide 
precipitation controlling strength. The pioneering work of Speich [30] established that good toughness in 
secondary hardening alloy steels demands near-completion of M 2 C alloy carbide precipitation in order to 
dissolve transient cementite (FesC) particles that otherwise contribute to microvoid nucleation during ductile 
fracture. The challenge to maintaining strength in such an overaged condition is to maintain a fine carbide 
particle size. A thorough electron microscopy study of M 2 C precipitation in the commercial AF1410 steel [31] 
was performed in the doctoral research of J. S. Montgomery (now at ARL/WMD) in conjunction with Small 
Angle Neutron Scattering (SANS) studies by Weertman and coworkers [32], demonstrating coherent 
precipitation behavior consistent with recent theory of precipitation at high supersaturations [33,34]. In this 
regime the fundamental scaling factor governing particle size (and thus alloy strength) is the initial critical 
nucleus size which scales inversely with the thermodynamic driving force for precipitation; the time scale of 
precipitation is governed by a diffusion time defined by the coarsening rate constant. To compute the 
thermodynamic driving force for coherent M 2 C precipitation, the doctoral research of K. C. King with 
Professors T. Mura and P. Voorhees addressed the elastic self energy of the coherent carbides [35], and also 
explored the elastic interactions with dislocations during heterogeneous nucleation [36]. Incorporating the 
composition dependence of the carbide lattice parameters, the THERMOCALC chemical thermodynamic 
database [22] was then modified to incorporate the additional elastic self energy contribution to the 
thermodynamics of coherent precipitation. To control the time scale of precipitation, theories of coarsening 
kinetics in multicomponent systems were developed by Kuehmann & Voorhees [37] and Umantsev [38] to 
define a multicomponent diffusional rate constant. 

Our recent work has tested model predictions in a series of 16Co-5Ni-0.24C steels with the Cr and Mo contents 
represented in Figure 3 displaying a contour plot of predicted coherent precipitation driving force. An X-ray 
diffraction study of the precipitation hardening behavior of these alloys was initiated by Visiting Scientist S. 
Endo of NKK Japan and correlated with model predictions by Research Scientist C. J. Kuehmann [4]. Measured 
peak hardness showed a strong correlation with the coherent precipitation driving force after taking into account 
the "paraequilibrium" condition of prior precipitation of cementite under carbon diffusion control. 

Toward improved precision of the M 2 C coherent thermodynamics, the doctoral research of C. Knepfler [39,40], 
with ceramist Prof K. Faber, addressed the synthesis and characterization of (Mo,Cr,V,Fe) 2 C carbides in bulk 
form, measuring the composition dependence of lattice parameters, thermal expansion and isotropic elastic 
moduli as inputs into a more precise treatment of the composition and temperature dependence of the carbide 
elastie self energy. Measured elastic constants were also applied to an analysis of the critical particle size for the 
shearA)ypass transition in precipitation strengthening of steels. 

The doctoral research of R-H. Liarng [41] with Prof Mura incorporated Knepfler's measurements in refined 
calculations of the carbide elastic self energy in an iron-based matrix, taking into account the differing elastic 
constants of matrix and particle, and the role of particle shape. Liarng also computed internal stresses on carbide 
shear systems relevant to coherency loss mechanisms, predicted solute distributions around coherent particles. 





and considered interparticle elastic interactions promoting autocatalytic nucleation. As the transformation 
eigenstrains for a coherent M 2 C carbide in an Fe matrix are quite high, these linear elastic calculations are 
regarded as upper bound estimates of self energies and local stresses. To rationalize observed coherent 
nucleation conditions and coherent particle compositions measured by atom-probe microanalysis, it has been 
necessary to apply a correction factor of about 1/3 to the computed self energies. Liamg performed lower bound 
self energy estimates by excluding the energy density contributions of highly strained regions likely to be 
strongly nonlinear, similar to the dislocation core cutoff approximation, supporting correction factors of this 
magnitude. More recently, preliminary nonlinear numerical calculations by Graduate Fellow L. Coates working 
with Prof Brian Moran offered further support for the reduced elastic energy. 



Figure 3. Contour plots of predicted coherent M 2 C preeipitation driving force for 16Co-5Ni-0.24C steels at 510 °C vs. Cr and Mo 
content. Open points denote compositions of experimental alloy series. 
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Figure 4. Evolution of microstrain distribution width (measured with Fourier column length of L=144A) with 510°C tempering of 
16Co5Ni4Cr4Mo0.24C steel. Dashed curve depicts Vickers hardness (right scale). 


In our most recent work, Visiting Scientist Y. Nagataki of NKK Corporation collaborated with the late Prof. J. 
B. Cohen [42] to perform precision X-ray line broadening studies to quantify both dislocation recovery during 
tempering and the evolution of matrix micro-strain associated with carbide coherency. The latter is represented 
as a function of tempering time in Figure 4 for the case of the highest driving force composition (4Cr4Mo) of 
Figure 3. The dashed curve representing Vickers hardness (right scale) shows that maximum coherency micro¬ 
strain occurs just before peak hardness. From corresponding SANS measurements of particle size, this identifies 
a critical particle size for coherency loss of 3nm. 
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Treating the high-supersaturation precipitating system as an evolving unstable equilibrium, Research FelloAV 
Carrie Campbell [43,44] integrated the diverse experimental data to correlate particle size (corrected for shape 
dependence of capillary energy) with thermodynamic driving force defined by measured phase fraction, as 
represented in Figure 5 for the example of the 4Mo alloy of the series of Figure 3. Slopes define interfacial 
energies for the coherent (below 3nm size) and incoherent states of the M 2 C carbides, supporting an elastic 
energy correction factor of 0.4 (as described earlier) to give the linear behavior shown for the coherent state. 
With these parameters, the thermodynamic database can predict the full trajectory of particle size, phase 
fraction and phase composition during precipitation. 



Figure 5. Correlation of measured effective particle size (corrected for shape dependence of capillarity) with reciprocal driving force, 
defining interfacial energy for coherent and incoherent M 2 C carbides in 16Co5Ni4Mo0.24C steel. 
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Research Assistant John Wise [45,46] incorporated the precipitation trajectory model into a strengthening 
model by calibrating the contributions of solution strengthening (ass), particle strengthening (op), martensitic 
interfacial strengthening (a„), and dislocation substructure strengthening (ao) as represented by the equations 
listed in Table 1 and calibrated from Small Angle Neutron Scattering (SANS) data in Figure 6. The figure show 
the calibration of the strengthening contributions of shearable and nonshearable (Orowan) obstacles to the 
SANS measurements. Using the measured particle size distributions, the analysis gives a critical particle size 
for the shear to bypass transition of 3nm, equivalent to the critical size for coherency loss. Figure 7 summarizes 
the evolution of predicted hardness with relative phase fraction (normalized to the final equilibrium fraction) for 
four of the alloys of Figure 3. The predictions are in excellent agreement with experiment, and show that higher 





driving force compositions achieve greater peak hardness by reaching the critical shear/bypass size at higher 
relative phase fractions. 



{GpR/2jvG^b}^‘^/{L-2R) (A'^) ]n(2R/rg)/(L-2R) (A"’) 

Figure 6. Strengthening model equations employed for secondary hardening martensitic steels. Lower figures show calibrations of 
particle strengthening contributions of shearable and unshearable (Orowan) obstacles from SANS data. 

Research Associate Dr. Gautam Ghosh [47] has performed an analytical electron microscopy study 
demonstrating the paraequilibrium state of cementite (formed without substitutional diffusion) prior to M 2 C 
precipitation, and has quantified the early stages of cementite composition change during M 2 C precipitation. 
The evolution of cementite composition is important to evolution of the matrix carbon potential, which exerts a 
strong influence on the driving force controlling M 2 C particle size. Toward tighter control of precipitation 
strengthening, a quantitative modeling of cementite/matrix substitutional solute partitioning using the DICTRA 
multicomponent diffusion code is needed. 
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Figure 7. Predicted evolution of hardness w ith relative carbide phase fraction in secondary hardening 16Co5Ni0.24C steels of varying 
Cr and Mo content. 

The evolution of M 2 C particle composition during coherent precipitation is very important to the calibration of 
elastic self energy models. A new opportunity for more rapid measurement of particle composition at the 
important 1-3 nm size scale is the establishment of a high mass resolution 3D atom probe at Northwestern based 
on a new design from Oxford University commercialized by Kindbrisk, Ltd. Testing of the new instrument at 
Oxford produced the reconstructed 3D images in Figure 8 for the 0.24C-16Co-5Ni-4Cr-4Mo steel in the peak 
hardness condition. The image highlights matrix Fe, Co and Ni and the carbide Cr, Mo and C. With discrete 
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crystal planes clearly evident, the 3 nm carbide corresponds to the critical size for coherency loss and 
shear/bypass transition. 
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Figure 8. 3D Atom Probe reconstruction of MaC carbide in 0.24C-16Co-5Ni-4Cr-4Mo steel solution treatment 1000 ®C, Ihr, 
tempered 500 °C, 2hr to peak hardness of Rc 60. 
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2. Toughening 


Early work under the SRG Program addressed toughening of UHS steels via (a) control of microvoid nucleating 
particle dispersions governing plastic shear localization [48-51], and (b) dispersed-phase transformation 
toughening by precipitation of optimal stability austenite particles [52,53]. A continuing effort under ONR 
support has defined requirements to resist intergranular fracture [54]. Guided by model predictions on 
microvoiding behavior, the doctoral research of C. J. Kuehmann under ARO support addressed optimization of 
the solution treatment of the commercial AermetlOO alloy to obtain a residual grain refining dispersion of fine 
Tie particles with improved resistance to microvoiding as attested by Jic toughness measurements [55]. 
Kuehmann then explored optimization of multistep tempering treatments to achieve dispersed-phase 
transformation toughening in this alloy. 

Our understanding of transformation toughening mechanisms benefited from a parallel DOE-sponsored effort 
studying fully austenitic alloys, experimentally demonstrating the role of phase stability and dilatation [2, 56- 
58] with numerical simulations defining the effect of transformation plasticity in enhanced resistance to plastic 
localization [59] and microvoiding processes [60]. ARO sponsored research has explored the application of this 
understanding to the control of dispersed-phase transformation plasticity in the secondary hardening martensitic 
steels. Our progress is summarized by the toughness-hardness plot in Figure 9. Early experiments guided by 
thermodynamic predictions of achievable austenite stability in the AF1410 composition boosted toughness to 
the AF1410-TT band employing multistep tempering to nucleate a fine dispersion of Ni enriched austenite 
while maintaining a sufficiently fine carbide dispersion for desired strength [52]. 
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Figure 9. Toughness-hardness plot with dark bands denotes properties of transformation toughened (TT) materials; dashed box 
denotes original SRG objectives. 

Our efforts to apply the same strategy to the higher strength AermetlOO alloy revealed a much tighter 
processing window due to the faster carbide precipitation kinetics associated with this higher Cr composition. 
Initial multistep tempering studies by Kuehmann [55] demonstrated transformation toughening, but at less than 
desired strength levels. A more thorough study was then undertaken by Visiting Scientist M. Srinivas from the 
Defence Metallurgical Research Laboratory, India, establishing first the detailed low temperature tempering 
kinetics for optimal strengthening, then interposing short intermediate temperature austenite nucleation 
treatments. Detailed measurement of toughness evolution was provided by relatively simple crack-tip stretch- 
zone measurements followed by verification of promising treatments with Jic tests. The results of optimal 
processing are represented by the AermetlOO-TT and MTLl-TT band in Figure 9. The outstanding toughness 
and hardness combinations are seen to lie within the dashed box representing the original property objectives of 
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the SRG program. Under a one-year contract from ARL/WMD (formerly MTL), SRG design tools were applied 
to the design of a prototype high hardness armor steel. This led to the alloy MTLl, which was evaluated by 
Kuehmann and J. Cho [55]. Employing a V addition to enhance strengthening efficiency, the alloy 
demonstrated properties equivalent to the AermetlOO alloy, but with lower C content (0.21C). Ballistic testing 
at ARL/WMD demonstrated V 50 levels at fixed plate thickness equivalent to the best monolithic armor steels. 

Austenite precipitation behavior in these alloys was then investigated by electron microscopy and microanalysis 
in the doctoral research of fellowship student H. Lippard [61] cosupervised by Prof. V. Dravid. Dark-field 
microscopy revealed two forms of precipitated austenite in the AermetlOO-TT material. The microstructure is 
similar to our previous observations in AF1410-TT [52] in which slightly Ni enriched austenite discs at 
martensite lath boundaries coexist with a fine dispersion of more Ni rich intralath precipitates; the lath boundary 
precipitation gives nonequilibrium austenite compositions while the composition of the fine intralath 
precipitates is in reasonable agreement with equilibrium calculations. The high austenite stability responsible 
for toughening these UHS steels is attributed to the fine size and Ni enrichment of the latter. While the fine 
austenite in AF1410-TT had a particle size of 20nm, the fine particles observed in AermetlOO-TT are only 5nm, 
taxing the limits of analytical electron microscopy. In addition to the 200 kV cold field emission Hitachi 
HF2000 analytical electron microscope at Northwestern, the samples were analyzed with a Phillips CM20 at 
Gatan Laboratories (employing a Gatan Imaging Filter for EELS imaging with the Ni core loss L-edge) and the 
highest brightness SOOkV cold field emission UHV VG603 STEM at Lehigh University. The high resolution 
microanalysis quantified the fine particle compositions throughout the precipitation process, and revealed 
evidence of austenite nucleation on the Cr-rich carbides. Consistent with the carbide capillary energies defined 
by the analysis of Figure 5, Lippard extended the unstable equilibrium model to treat the trajectory of co¬ 
precipitation of austenite and M 2 C. Figure 10 presents computed austenite Ni contents as a function of capillary 
energy in both austenite and M 2 C. The boxes depict the trajectory of co-precipitation indicated by the austenite 
composition measurements. A greater than 50% drop in the M 2 G capillary energy prior to austenite precipitation 
is consistent with austenite nucleation on the carbides. 

While the analytical electron microscopy gives a plausible representation of co-precipitation, the new 3D atom- 
probe capability offers the opportunity for definitive measurement of both the austenite and associated carbides 
during austenite nucleation. A thorough investigation with this technique and an extension of the carbide 
precipitation computational model will allow full simulation of the co-precipitation of austenite and M 2 C. 
Although the toughening results represented in Figure 9 demonstrate the potential of transformation toughening 
in these steels, the optimal austenite nucleation treatments in both AF1410-TT and AermetlOO-TT required 
tempering times of 10 to 15 minutes, which would not be practical in many applications. A reliable, predictive 
co-precipitation model will be essential for the robust design of alloys capable of reproducibly achieving 
optimal transformation toughening with practical heat treatment conditions. 

3. Gear Steels 

Case hardenable secondary hardening steels for the high power density gear systems of Figure 2 may represent 
the most promising area for application of the design models we are developing. This line of research originated 
as conceptual designs in Materials Design class projects performed in collaboration with a mechanical 
engineering design class who identified property objectives for a 50% weight reduction in gears. Preliminary 
designs comparing low temperature nitriding and high temperature carburizing led to selection of carburizing as 
the most promising processing route to achieve the Imm-scale case depths needed in typical gear applications. 
Two generations of carburizable prototype compositions were designed and evaluated. An undergraduate 
project evaluating the second prototype alloy, designated C2, won 1st prize in the TMS-AIME national student 
design competition [ 21 ]. 
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M^C Capillary Energy (J/mol) 

Figure 10. Contour plot of precipitated austenite Ni content (wt. pet.) for AermetlOO steel tempered at 455 °C, showing observed 
trajectory of austenite precipitation vs, M 2 C and austenite capillary energy. 



Figure 11. Probability distributions of RCF fatigue life in standard NTN ball-on-rod RCF tests comparing ultrahard case steels with 
M50 bearing steel. 
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Research on carburizable gradient systems for gear applications continued in the doctoral research of J. 
Wise[45,46], initiated under seed funding from the Ford Motor Company. Precise treatment of the dynamics of 
gradient formation by carburizing has been made possible by the DICTRA multicomponent diffusion code. To 
apply the code to carburization of the complex multicomponent steels of interest, Wise conducted preliminary 
carburizing experiments on an Fe20Col0Ni alloy and determined the composition-dependent C diffusivity from 
measured C profiles. Based on these measurements, it was found necessary to modify the DICTRA mobility 
database (mob4) through a Co-C interaction term. The model was then validated by comparison with measured 
diffusion profiles in carburized AermetlOO steel, demonstrating excellent agreement. Combining DICTRA 
diffusion simulations with the strengthening model of Figure 7 provides the basis for precise control of hardness 
profiles in these secondary hardening steels. The hardness profiles so far obtained in prototype gear steels 
compared with that of a conventional gear steel, demonstrate a case hardness increase of 50%. The 
corresponding increase in rolling contact fatigue (RCF) life obtained thus far is shown in Figure 11, based on 
the standard NTN ball-on-rod test. When carburized to case carbon content above the design level of 0,70C, 
fatigue life is comparable to the M50 bearing steel, and metallographic examination indicates fatigue initiation 
at large undissolved cementite particles. When the case is properly diffusion treated to the desired 0.70C level, 
these excess carbides are eliminated, and preliminary RCF data (with and without CrN coatings) indicates a 
narrower probability distribution with an order of magnitude life improvement over M50 at a 10% failure 
probability. 

Research Assistant Jim Wright is continuing the RCF testing with a higher cleanliness heat of the experimental 
steels, and is investigating the nature of the microstructural instabilities causing "butterfly" and "white band" 
defects that precede RCF initiation in conventional gear and bearing steels. Wright is also performing X-ray 
diffraction studies to investigate the control of case residual stress in the new secondary hardening steels. 

New developments in numerical applied mechanics methods provide a unique opportunity for fundamental 
modeling of the multiscale phenomena underlying the microstructural instabilities driving RCF initiation The 
multiple structural scales of concern are represented schematically in Figure 12. A typical contact width for a 
bearing is about 0.1 mm and contact stresses are significant over a depth (h) of about 1mm, The stresses 
induced by contact interact with subsurface inclusions which are about 2 microns in diameter. The material 
behavior on the submicron scale in the region of local stress concentrations associated with these inclusions is 
governed by the mechanisms of kinematic hardening and microvoid softening associated with 0 . 1 pm scale 
second phase particles such as M 3 C carbides. 



Figure 12. Multiple structural scales in RCF initiation represented for rolling element bearing. 
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The kinematic hardening is due to the Bauschinger effect associated with the back-stress induced by the second 
phase particles. This kind of behavior can be modeled using kinematic hardening plasticity models. Aimer, 
Cohen and Moran [62] carried out finite element simulations of cyclic plasticity using a generalized plasticity 
model with kinematic hardening. The kinematic hardening was due to residual microstresses in ferrite induced 
by the inhomogeneous deformation fields around cementite particles. The finite element method was used to 
track the local stress-strain behavior in the vicinity of a stress concentration during a variety of press-fit 
operations. Aimer et al. found good agreement between calculated residual stresses and those obtained by X-ray 
measurement techniques for dual phase materials. An interesting feature of this work is that the kinematic 
hardening parameter, for use in the constitutive model, was determined from X-ray measurements. 

The microvoid softening is due to the nucleation and growth of microvoids at the second phase particles, which 
can ultimately lead to void-coalescence and the formation regions of strain localization or shear bands. This 
kind of behavior can be modeled by a Gurson type constitutive relation, which incorporates the nucleation and 
growth of voids. 

4. Biomimetlc Laminates 

Testing the limits of complexity in rational materials design Incorporating biomimetic adaptive system 
concepts, the essential concepts of the self healing "smart steel" superalloy composite are illustrated in Figure 
13 and the flow block diagram, which summarized the concepts as a system chart, is shown in Figure 14. A 
ferritic (BCC) superalloy strengthened by a coherent dispersion of ordered B2 precipitates is reinforced by a 
thermodynamically compatible y strengthened austenitic (FCC) shape memory alloy in a laminate formed by 
hot rolling of a consolidated blend of prealloyed powders. The composite integrates two biomimetic concepts. 
At low temperatures the shape memory alloy reinforcement toughens the otherwise brittle ferritic alloy by 
bridging cracks via martensitic pseudoelasticity, mimicking the nonlinear polymer reinforcement of a brittle 
ceramic in the toughened laminate system of natural seashell [10,11]. Returning an internally damaged 
composite to the superalloy operating temperature, the bridges contract via the shape memory effect providing a 
crack clamping force to promote diffusional rewelding of cracks as a biomimetic self-healing mechanism. 
Thermodynamic compatibility demands that the BCC (a) and FCC (y) alloys be in near-equilibrium with each 
other at hot working and solution treatment temperatures, and that the two 2-phase a+B2 and y+r' alloys reside 
in an equilibrium 4-phase a+B2+Y+Y' field at service temperatures. Matrix/reinforcement interfacial strength 
must be limited to promote crack bridging behavior, and the austenitic alloy must exhibit higher creep strength 
than the ferritic alloy to maintain crack clamping forces for damage healing. 




Figure 13. Schematic of biomimetic laminates. At low temperature, crack are formed, causing a thermoelastic martensite 
transformation in the reinforcement. Upon heating, the reverse martensitic transformation occurs, causing crack closure. 
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The original conceptual design and preliminary feasibility calculations by an undergraduate Materials Design 
class project team won 2nd prize in the first TMS-AIME national student design competition [20], Under 
AASERT support, the doctoral research of Brad Files [63-66] explored the design of a precipitation- 
strengthened NiAl intermetallic reinforced by a thermodynamically compatible y'-strengthened austenitic shape 
memory alloy. Thermodynamic models of the Fe-Ni-Al-Ti-Co-Mo system was refined with explicit treatment 
of order/disorder transitions underlying miscibility gaps. The models identified a suitable 4-phase a-B2-Y-Y' 
field at a service temperature of 600 °C in the basic Fe-Ni-Al-Ti system. Diffiision couples of alloys A and B, 
corresponding to a+B2 and y+Y* alloys with dispersion fractions of 20% were prepared and reacted 24 hrs at 
llOOC to investigate the diffusional interaction at solution temperatures where each alloy should be single¬ 
phase. The electron microanalysis data (employing a Hitachi S-4500 FEG-SEM) demonstrated a sufficiently 
large diffiision distance of 50-100|im allowing accurate measurement of multicomponent diffusion composition 
trajectories. The measured trajectories were consistent with the predicted high A1 diffusivity. The boundaries of 
the single phase fields at the respective Ti and A1 contents of alloys A and B were calculated based on our 
database. These calculations indicated that alloy A was actually well within the two-phase a+y field at 1 lOOC. 
To simulate a couple between two single-phase alloys, DICTRA simulations were run modifying the alloy A 
composition to the nearby single-phase a composition. TEM microscopy revealed a martensitic structure in the 
alloy indicating the alloy was fully austenitic at 1100 °C. These preliminary couple experiments provided 
important input for refining the thermodynamic database to properly position the relatively narrow y+a field in 
the region of interest, while also providing useful diffusivity information for further refinement of the DICTRA 
mobility database. The designed composite was fabricated as a laminate, demonstrating thermodynamic 
compatibility, crack-bridge toughening, and partial crack closure. 

In a feasibility study, TiNi shape memory alloy wires were used as the reinforcement in a Sn-Bi matrix to 
attempt to show that a healing behavior can be mimicked in engineered structure materials. When the composite 
was pulled to 1% strain, 80% of the plastic deformation was returned by heating above the reversion 
temperature of the shape memory alloy fibers. In the case of a composite strained above 6%, surface 
macrocrack growth was stabilized by fiber reinforcement. Upon heating this sample close to the matrix melting 
point, fiber reversion provided crack closure and clamping, with an accompanying reversal of 2% strain [64]. A 
Freshman Design team (Terminator 3) then designed a Sn-Bi-In alloy matrix with controlled melting 
characteristics. As a Sophomore project, Priscilla Bemikowicz fabricated the designed composite, demonstrated 
the predicted crack healing by exploiting incipient melting, and won the 1998 TMS-AIME undergraduate paper 
competition [67]. A ternary Sn-based matrix alloy reinforced by TiNi shape memory alloy was shown to be 
capable of closing cracks and healing them when heated to a partially liquid state [68]. The roles of porosity, 
composite wire content, and phase fraction liquid during healing were assessed using metallography as well as 
the Thermo-Calc thermochemical software system. It was shown that a strength recovery of 88% was achieved 
in the cracked samples after healing. Following the TMS-AIME award winning work, the composition of the 
matrix alloy was further refined to increase the temperature range over which the alloy can be heated and a flux 
was designed to reduce oxidation of the sample and to increase the flow of the liquid alloy [69,70]. 

Over the past few years, several Terminator 3 prototypes have been developed by various student teams from 
the Engineering Design and Communications (EDC) program as well as from the MSc 390 Materials Design 
classes at Northwestern University. It has been demonstrated that this project is not toO far away from bringing 
fiction to life. The concept of self-healing and “smart materials” has drawn attention of NASA-Houston, 
interested in having materials that can self-heal for remote repairing of space station parts. To make the alloy 
even smarter, a control system, or “brain” was designed to detect cracks and to command the heating procedure 
[71]. To accelerate the healing process the shape memory behavior of the TiNi alloy (Nitinol) was further 
investigated [72]. It was found that aging plays an important role in Ti-rich Nitinol wires. T^Ni phase is usually 
found in the alloy, lowering the shape memory temperatures. Solution treatment can dissolve significant amount 


of Ti 2 Ni particles. Increasing the Ti content, the Af temperature of the alloy (martensite to austenite 
transformation temperature) increases about 30 °C. The healing process is thus accelerated at the higher healing 
temperature. A complete crack healing was observed within 3 hours. 
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Figure 14. System structure of self-healing "smart steel" superalloy composite. 


Under continuing AASERT support, doctoral student Deborah Burton, working with Prof Moran, has 
implemented shape memory constitutive laws in a finite element simulation of such shape memory reinforced 
composites. The model will be applied to the prediction of optimal matrix and reinforcement characteristics for 
the extension of the healing concept to higher strength materials for application in smart structures. 
















































Conceptual Designs: Metal/Ceramic Hybrid Armor 



A conceptual design activity led by Dr. Jian-Sheng Wang addressed various approaches to hybrid metal ceramic 
armor materials design beginning with the system block diagram of Figure 15 representing possible functions of 
the primary structural subsystems. One possible structural configuration considered in greatest detail is 
represented in Figure 16. In accordance with the requirements of figure 15, the system consists of 5 layers. The 
front layer is an organic coating for the protection of environment attact, confinement of ejecta and possibly 
others such as antiradar function etc. Our research has not assessed this layer in detail. The second layer is a 
hard layer for projectile breakup, deflection, erosion and load spreading. The requirements for this layer are: 
ultra high hardness combining with high toughness, constrained fracture and ability of crack arrest to limit 
impact damage and formation of debris, and capable of spreading impact load to reduce ballistic penetration. 
We propose a transformation-toughened UHS alloy grid with bonded CMC ply tiles as the hard layer. In order 
to reduce the vehicle weight, we explored the feasibility of transformation-toughened Ti-based alloys. The third 
layer is a possible nonlinear compliant material (rubberlike) for the damping of vehicle dynamics and 
absorption of ballistic impact energy. The fourth layer is a smart ceramic layer consisting of an SPPS ceramic 
matrix with Sn-based low melting-point particles and SMA reinforcement fibers. The SPPS ceramic matrix 
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provides ballistic penetration resistance. The SMA fibers provide bridging toughening under impact loading and 
crack closure and clamping effect during self-repairing. The Sn-based low melting-point particles provide self- 
healing for cracks. An acoustic sensor is attached to the layer to detect the damage. Whenever damage is 
detected, a high density low voltage DC current will be automatically activated through the SMA fibers heating 
up the Sn-based low melting-point particles providing self-healing of cracks. The fifth layer is a load bearing 
layer made of a case-hardened functional gradient UHS steel. This layer is also the substrate of the SPPS thick 
coating. The hardness varies gradually from the surface to the core. The high hardness surface provides further 
ballistic penetration resistance, the tough core provides load bearing, energy absorption and debris confinement. 
To restrict crack propagation, a laser local annealing technique could be applied to form a patterned soft 
network providing crack arrest function. 
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Figure 16. Schematic structure of an integrated multifunctional hybrid armor system. 

1. The CMC ply tiles 

An ultra high hardness is the primary requirement against ballistic penetration, this hard layer has to be ductile 
to resist debris formation. A single phased ceramic tile does not satisfy this requirement. The central technical 
problem for “toughened” materials is their lack of ductility and the associated scale dependence caused by 
weakest link dominance and the extreme value nature of the stochastic behavior. The weakest link scaling can 
be suppressed by imparting damage tolerance to the material, through mechanisms that stabilize sub-critical 
cracks. For this purpose, we propose CMC cross-plies as the anti-ballistic tiles. The advantage of the CMC 
cross-ply is that it provides toughness by the mechanisms of fiber pullout and crack deflection while 
maintaining the high hardness of the ceramic. We propose a 90 degree fiber reinforced CMC ply. The 
toughening effect induced by inelastic strains has been documented [73-76] and schematically shown in Figure 
17 [77]. Upon tensile loading, cracks form in the 90° plies by tunneling with associated inelastic strain (Figure 
17). Upon subsequent loading, as the cracks penetrate the 0° plies, they interact with the fibers and the coating. 
When debonding and slip occurs within the coating, two stresses characterize the inelastic strain: a friction 
stress T, and a debond stress Gi (Figure 17). The latter is related to the debond toughness for the coating and the 
residual stress. The consequent inelastic tensile strain depends on the stress acting on the 0° plies. 
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We propose TiB 2 as the matrix of the CMC plies. A comparison 
of thermomechanical properties of various ceramics, including 
oxides, carbides and borides, is shown in Table 2. TiB 2 shows the 
best combination of the ultra high hardness with good toughness. 

The hardness of TiB 2 is less than diamond and cubic BN, ranking 
as the third hardest material in nature. Considering the primary 
candidate materials of the confinement grid, Ti-based alloys, 
there are other advantages of choosing TiB 2 as the matrix 
material: the crystallographic structure, the density, and the 
coefficient of thermal expansion (CTE) of TiB 2 are comparable 
with Ti-based alloys. It is expected that the residual stress at the 
CMC/Ti interface is low, and the bonding between them is 
strong, inducing large confinement due to crack arrest and load 
spreading. 

Figure 17. The ply cracking model and the associated stress-strain curve (a) 
and the cell model for inelastic strains in the 0° plies. 

Commercially available ceramic fibers are limited to Nicalon and 
Tyralon, which are SiC-based fibers. SiC-matrix CMCs with 
Nicalon or Tyralon reinforcement fibers and BN coatings has 
been extensively studied. We propose use of these commercial fibers in our TiB 2 matrix CMC. We next need to 
investigate the compatibility between the fiber and the matrix, the thermomechanical behaviors of the CMC, the 
processing and the anti-ballistic performance of the CMC tile. 

2. Constraint Frame 

We propose transformation-toughened UHS steels and Ti-base alloys as the constraint grid materials to provide 
constraint for crack propagation in the CMC tiles and preventing debris. Our imderstanding of dispersed-phase 
transformation toughening mechanisms has been summarized earlier. Plastic deformation causes retained 
austenite particles to transform to martensite inducing a delayed hardening effect and enhanced resistance to 
plastic localization. ARO sponsored research has explored the application of this understanding to the control of 
dispersed-phase transformation plasticity in the secondary hardening martensitic steels. Outstanding toughness 
and hardness combinations have been achieved in AermetlOO-TT and MTLl-TT steels. 

Table 2: _ 
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2.8 
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The same principles for UHS steels can be applied for the investigation of phase-transformation toughening Ti- 
based alloys. The choice of Ti-based alloy as the constraint grid material can reduce the weight of the vehicle 
and improve the compatibility between the grid material and the CMC tiles. We propose design of a phase- 
transformation toughening UHS Ti-based alloy. Ti-based alloys can undergo martensitic transformation the 
same as in steels. THERMOCALC/DICTRA modeling can explore the feasibility of developing a Ti-based 
alloy with retained “austenite” particles,, addressing precipitation hardening, the stability of the retained 
austenite and the kinetics of the stress induced martensite transformation in the alloy. 

3. The smart material 

We propose a smart ceramic composite thick coating consisting of an SPPS ceramic matrix with shape memory 
alloy (SMA) reinforcement fibers. The composite integrates two biomimetic concepts: at room temperature the 
shape memory alloy reinforcement toughens the ceramic by bridging cracks via martensitic pseudoelasticity, 
mimicking the nonlinear polymer reinforcement of a brittle ceramic in the toughened laminate system of natural 
seashell. At an elevated temperature, the bridges contract via the shape memory effect providing a crack 
clamping force. To promote diffusional rewelding of cracks as a biomimetic self-healing mechanism low 
melting-point Sn-based alloy particles can be co-deposited by SPPS processing. To make the ceramic even 
smarter, a sensor for detecting cracks and an actuator for initiating heating of the ceramic can be attached. 
Whenever damage of the system is detected, a high density, low voltage DC current can be automatically 
activated through the SMA fibers heating up the Sn-based low melting-point particles providing self-healing of 
cracks. 

The matrix of the smart ceramic composite we proposed is partially stabilized Zr02 (YSZ). YSZ thick coating 
can co-deposit with a self-healing Sn-based alloy by SPPS process. Reinforcement fibers of TiNi SMA (Nitinol) 
can be added during deposition. YSZ has the highest toughness among the ceramics listed in Table 2. The 
comparable elastic modulus of YSZ with that of SMA insures a easy crack closure by contraction of the SMA 
fibers. The extremely low thermal conductivity of YSZ insures efficient heating up of the smart ceramic for 
self-healing. The thermodynamic compatibility between the components must be addressed along with the 
thermomechanical properties, the bridging toughening effect and the self-healing behavior of the smart ceramic. 
With an acoustic device attached to the system, one can study the sensibility of the acoustic emission in 
determining the crack density and determine the critical value of the acoustic signal under the pre-determined 
damage tolerance. The critical signal can be sent out to actuate a switch to allow a low voltage-high density 
current passing through SMA fibers heating up the composite to the Ai temperature of the SMA. At this 
temperature, reverse phase transformation occurs in the deformed SMA and the SMA fibers return to the 
undeformed shape inducing crack closure of the damaged composite. To insure rewelding of the crack surfaces, 
particles of a low melting point alloy are co-deposited with YSZ by SPPS processing. The melting point of the 
alloy should be slightly lower than the Ai temperature of the SMA. The alloy must have a good wetability with 
YSZ and a low viscosity at the liquid state. The Terminator 3 project can continue to design the Sn-Bi-In alloy. 
The composition of the alloy, the roles of porosity and fiber content of the composite need be assessed using 
metallography and the THERMOCALC/DICTRA software system. 

4. Strengthening of a functionally graded case hardened steel 

We propose a case-hardened functional gradient UHS steel as the load bearing layer of the armor system. The 
hard surface provides further ballistic penetration resistance, the tough core provides load bearing, energy 
absorption and debris confinement. To restrict crack propagation in the hard layer, a laser local annealing 
technique can be applied to form a patterned soft network providing crack arrest function. This load bearing 
layer is coated with the smart ceramic coating. 
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High temperature carburizing has been shown the most promising processing route to achieve about 1 mm-scale 
case depths needed. Employing a V addition to enhance strengthening efficiency, the MTLl alloy has 
demonstrated properties equivalent to the AermetlOO alloy, but with lower C content (0.21C). Ballistic testing 
at ARL/MD demonstrated V 50 levels at fixed plate thickness equivalent to the best monolithic armor steels. A 
higher carbon composition (0.25C) may further increase the hardness. This prototype alloy may be too costly, 
however, for applications in the armor system. A new alloy with less content of Co and Ni can be designed by 
THERMOCALC computations. Co could be replaced by Al, and Ni can be partially replaced by Mn. The 
THERMOCALC/ DICTRA code can be used to design carburization of the new complex multicomponent 
steels. Combining DICTRA diffusion simulations with the strengthening model provides the basis for precise 
control of hardness profiles in these secondary hardening steels. 


Exploration of Shearable Matrix Alloys for Tungsten Heavy Alloy Penetrators 

Based on the premise that the shearability of U penetrator alloys is based on the transformation plasticity 
associated with adiabatic transformation to the softer bcc phase, the Senior Project of Andrew Elliott, aided by 
thermodynamic modeling of Dr. Gautam Ghosh, addressed design of Hf-Ta-based alloys with identical 
transformation characteristics, for application as matrix alloys for shearable WHA penetrators. After 
thermodynamic assessment of the Hf-Ta binary system, button alloys were prepared by both elemental powder 
sintering and remelting. Although transformation temperatures could not be clearly measured, alloy hardness 
levels of 330 VHN were obtained with desired martensitic microstructures. Tests at ARL have shown that Hf- 
Ta alloys demonstrate desired adiabatic shear behavior. A marked copy of Elliott’s report is attached as 
appendix B. 
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Appendix A 


PATHWAYS OF DISCOVERY 

Designing a New 
Material World 

Gregory B. Olson 


Materials have paced the evolution of technology for millennia. Their importance in the advance of 
human civilization is apparent in the naming of historical epochs, from the Stone Age through the 
Bronze and Iron Ages and into the ongoing Silicon Age. The origin of diversity in the material 
world remains largely mysterious to the public, yet the specialists’ ability to understand and manip¬ 
ulate the microstructures of materials has grown explosively over the past half-century. As the new 
millennium unfolds, a confluence of natural philosophies—one that combines reductionist and syn¬ 
thetic thinking—is ushering in an Age of Design marked by new materials and ways of creating 

them that go beyond the dreams 
of the medieval alchemists. 

Materials as Systems 

The modern view of materials 
structure was best expressed by 
the late philosopher-scientist Cyril 
Stanley Smith ( 1 ). He described a 
universal multilevel nature of 
structure with strong interactions 
among levels and an inevitable interplay of perfection 
and imperfection at all levels. Smith argued that the 
materials scientist’s distinct view of structure is de¬ 
fined by the desire to understand the structure and 
property relations underlying the technological and 
economic value of materials. 

This view of matter integrates science and engi¬ 
neering and is built on a natural philosophy that is 
older than science. Smith identifies the origin of this 
philosophy within medieval alchemy, whose practi¬ 
tioners are mostly remembered for their attempts to 
“transmute” base metals into gold. They knew noth¬ 
ing about protons and electrons. Still, Smith asserts that their aspirations, motivations, and even 
concepts resembled the modem materials engineer’s use of complex processing to “transmute” the 
multilevel microstmcture of materials (from the atomic to macroscopic levels) to achieve the es¬ 
sential property of gold, namely, economic value. 

The property-driven view of structure and processing for the creation of value, shared by modem 
materials science and alchemy, embraces an essential complexity of material structure. During the 
17th century birth of modem science. Smith identifies a short-lived Golden Age of materials science 
under the leadership of Rene Descartes, whose “corpuscular” philosophy inspired the development of 
a multilevel view of stmcture to account for the diverse properties of materials. 

The prescient grasp of materials achieved by the corpuscularians is well represented by Rene An¬ 
toine Ferchault de Reaumur’s 1722 sketch of the stmcture of quench-hardened steel (see figure 
above). He proposed that a single grain of steel (G), if enlarged, would reveal a set of “molecules” 
(M) and voids (v). Higher magnification would reveal a substmcture of the molecules (pp); and yet 
higher magnification would show an aperiodic arrangement of spheres. The finest scale Reaumur en¬ 
visioned corresponds to a periodic pacing of spheres, what we might think of as the nanoscale. 

This elegant view contained many stmctural elements of modem materials science. There was no in¬ 
strumentation with which to validate it, however, so the complex stmctural view was supplanted by more 
intellectually compelling but overly simplistic notions. One was Isaac Newton’s continuum, which 
erased stmctural considerations entirely. The other was John Dalton’s atomism, which held that there 
was only one level of stmcture that mattered. These conceptual idealizations were sufficiently com¬ 
pelling to divert the corpuscularian framework for 2 centuries. 



Materials savant. Rene de Reaumur and his 
1722 sketch of steel’s anatomy. 
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We have had the past century to reinvent materials sci¬ 
ence. The atom and the continuum remain dominant philo¬ 
sophical forces; they are the foundation of existing theoreti¬ 
cal tools. But armed with the growing reser¬ 
voir of structural facts gleaned from instru¬ 
mentation and computation, the challenge 
now is to adapt these theoretical tools to un-. 
derstand and control the complex structures 
of real materials. 

The modem view of material stmcture 
differs from Reaumur’s mainly in the de¬ 
tailed morphologies characteristic at the 
different length scales of a material’s hier¬ 
archical stmcture. In the case of steel, the 
most significant difference is the over¬ 
estimation of porosity in the 18th century 
depictions. And if Reaumur’s voids are reinterpreted as 
free volume,” his sketch becomes a reasonable model of 
polymeric materials. It shows the remarkable ability of the 
human mind to infer necessary stmcture from the contem¬ 
plation of properties alone. 

Smith made another important historical observation 
about materials de¬ 
velopment (2). Since 
prehistory, people 
have put newly dis¬ 
covered materials to 
practical use long be¬ 
fore they understood 
much about them. 

Consider the pat- _ 

tern-welded sword edge. Etched and repolished tip 

blades made and 3 century sword, 
used by Merovingian 

Franks and Vikings (see figure above) as early as the 6th 
century A.D. By hammer-welding steels of differing car¬ 
bon content (a technical feature unknown to the early 
sword smiths), a laminate composite was fabricated by 
labor-intensive “hand lay-up” to produce a hybrid stmcture 
with a tough core supporting a hard cutting edge. 

The actual mechanism of the quench hardening of steel, 
which is responsible for the hardened sword edge, was hotly 
debated even in the early 20th century. There were two fac¬ 
tions: the “allotropists,” who favored stmctural transforma¬ 
tion, and the “carbonists,” who pegged the hardening mech¬ 
anism on dissolved carbon. A cmcial observation was the 
discovery by Floris Osmond in 1893 of “martensite” in 
steel, the microstmctural form that iron assumes during the 
quenching process. We now know that the combination of a 
martensitic transformation and the redistribution of trapped 
interstitial carbon (along with ambient aging during which 
additional microstmctural evolution occurs) underlies the 
edge hardness of the 6th century swords. 

The ancient swords embody another illustration of how 
art and craft has traditionally preceded science. The sword 
smiths’ use of chemical etching to bring out aesthetic metal¬ 
lic patterns set the foundation for modem metallographic 
observation of microstmcture established by Henry Sorby in 
the 19th century. By etching metal samples with acid, Sorby 
revealed internal microstmctures and correlated them with 
the properties and performance of the materials. 

Transmission electron microscopy (TEM), a modem ana¬ 
lytical cohort of metallography that reveals finer stmctural 
levels (see figure at right), suggests that the. ancient sword- 
makers ever were accidental nanotechnologists. TEM analy¬ 


sis reveals nanometer-scale patterns of carbon in the hard 
edge of the sword. These patterns emerged from a process 
called “spinodal decomposition”, in which a solid solution, 
such as the sword’s high-carbon martensitic steel, becomes 
unstable and its constituents reorganize (3). In modem par¬ 
lance, the swordmakers’ products could be described as “self- 
assembled heterophase nanostmctures.” 

The spatial dimension of materials is only part of what 
makes them tick. There is a spectrum of characteristic relax¬ 
ation times associated with the various chemical and physi¬ 
cal processes operating at the material’s differing stmctural 
length scales. This adds the dimension of time. 

The resulting dynamic spatiotemporal hierarchy means 
that any material at any time has stmctural features, such as 
grain sizes and dispersed particles, that have not yet 
reached equilibrium. That is why a material’s stmcture and 
properties depend on how it was made and what conditions 
it endured in service. 

Also fundamental to this dynamic view of materials is a 
realization that stmctural defects play vital roles at all length 
scales. Defects can make things fail in a thousand ways, but 
they’re also often what makes materials so valuable. Perfect 
_ silicon crystal without dopant ions is not the semi¬ 
conductor that has changed society. Pure iron met¬ 
al without the right spicing of carbon would never 
have become the steel backbone of the industrial 
revolution. Although “defect tolerance” remains a 
central tenet of modem materials science and is of 
incalculable commercial and safety importance, 
“defect engineering” is ascendant in the minds of 

_ many materials researchers. That’s because defects 

region of various hierarchical levels are a principal op¬ 
portunity for controlling material behavior. 

Because the personality of each material de¬ 
pends on all of these interacting spatial and dynamic at¬ 
tributes, it makes sense to approach materials as complex 
systems. Smith advocated such a systems view of materi¬ 
als decades ago. As more contemporary practitioners live 
by that insight, they are finding pathways to important 
new materials that can catalyze advances in manufactur¬ 
ing tools, computers, communications systems, and the 

_ myriad technologies 

whose very existence 
or improvement de¬ 
pends on more capa¬ 
ble materials. 




The Materials Disci¬ 
pline Comes of Age 

Two principal branch¬ 
es of natural philoso¬ 
phy have evolved to 
form modem materials 
Inside steel. TEM reveals nanoscale science. One is reduc- 
structure of a quench-hardened steel, tjonigt analysis, which 

takes nature apart to 
discern and understand its fundamental units. Reductionism 
has operated throughout the development of science. More 
sporadic has been the evolution of the synthetic systems 
view, which is better suited for understanding the connec¬ 
tions holding mature together. In a new balance of these two 
philosophies, the systems view integrates the fruits of our 
investment in reductionism while replacing conventional 
discovery-based R&D with a far more effective and efficient 
design-based approach. 


994 


12 MAY 2000 VOL 288 SCIENCE www.sciencemag.org 


CREDITS: (TOP) DE RE METALUCA; (MIDDLE AND BOTTOM) C. B. OLSON 







CREDITS: (TOP TO BOTTOM) C. B. OLSON; CO RBI S 



Planning materials. Flow-block diagrams guide materials design. 


Materials by Design 

To design a material is to try to meet a material 
user's need. A good place to start Is with proper¬ 
ty cross-plots, like those of Michael Ashby, that 
graphically define property-performance rela¬ 
tions (12). These help engineers select materials 
for their product designs. They're useful for 
defining a quantitative set of property objectives 
that will sum into the materials performance 
needed by users. These performance specifica¬ 
tions are determined by the role of the potential 
material in the wider system it serves. They also 
help define economic parameters, such as the 
cost of raw materials and processing, for the 
overall material design task. 

With objectives and economic constraints de¬ 
fined, the linear three-link framework of materi¬ 
als science and engineering (see figure, p, 996) serves to 
guide the design and development phases.The Steel Research 
Group (SRG) (see main text) uses flow-block diagrams that 
represent 1) key microstructural subsystems (such as ciystal 
grain sizes), 2) the primary links of these subsystems to the 
properties they control (such as strength and toughness), and 
3) the stages of processing (such as tempering or reheating) 
that govern their dynamic evolution. 

With these in hand, systems analysis is then applied to 
identify and prioritize key structure-property and processing- 
structure relations. Often, part of this exercise involves some 


additional modeling or empirical data gathering to fill gaps 
in the knowledge required for making practical decisions 
about composition and processing details. The systems view 
operates here at the strategic level, but It is supported by 
traditional reductionist analysis at a tactical level. On bal¬ 
ance, the procedure greatly reduces the amount of costly 
experimentation in materials creation. Instead of making 
tens of prototypes along the way to a useful new material, 
SRG designers reach their target metal using only a few ac¬ 
tual melts to refine the computation-heavy design efforts. 
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Although there were seeds of this turning point in materi¬ 
als science in early industrial laboratories at General Elec¬ 
tric, Bell Labs, and elsewhere, the multifaceted field of ma¬ 
terials research was deliberately synthesized as a single aca¬ 
demic discipline in the late 1950s and l'960s with the foimd- 
ing of the first university materials departments (4). Special¬ 
ists in the science and technology of metals, ceramics, poly¬ 
mers, and composites collaborated in pursuit of unifying 
principles for the creation of materials of all classes. This 
meeting of minds helped lay the sociological and cognitive 
groundwork for a systems approach to materials. 

The development of operations research in World War n 
and large-scale national missions such as the Manhattan Pro¬ 
ject and the Apollo space program were also important. It is 
estimated that at least 70% of our unprecedented economic 
boom of the past decade derives from technology, which in 
turn derives in good part from radically improved productiv¬ 
ity via new systems-based methods of product development. 

After World War H, newly formed government agencies, 
including the Office of Naval Research, the Defense Ad¬ 
vanced Research Projects Agency, the National Science Foun¬ 
dation (NSF), and o^ers also helped midwife the multidisci¬ 
plinary sensibihty needed for developing complex materials. 
The same was true in academe. Starting with Smith’s Institute 
for the Study of Metals at the University of Chicago, a series 
of university materials research laboratories estabhshed a na¬ 
tional infrastructure for the interdisciplinary enhancement of 
materials science. The private sector added sociological foun¬ 
dations for the emerging discipline through a diverse mix of 
professional materials societies, as well as the estabhshment 
of numerous materials-centered journals and conferences. 

The field now consists of many thousands of practition¬ 
ers who share a growing sense of community, yet whose col¬ 


lective expertise runs a stunning gamut of materials cate¬ 
gories. There are synthetic diamond makers, metal alloy de¬ 
signers, polymer scientists, optical fiber experts, thin flhn 
ceramic makers, developers of compound semiconductors, 
and “biomimetic” materials researchers who aim to emulate 
or adapt biology’s unmatched brilliance in materials innova¬ 
tion. The list goes on and on, with a diversity akin to the liv¬ 
ing kingdom’s millions of species. 

Over the past 4 decades in materials R&D, however, there 
has been a consistent emphasis on “good science,” as defined 
by reductionism, rather than “good materials,” which emerge 
when engineering, manufacturing, and economic factors are 
included in the mix. This bias has limited the technological 
impact of the materials science community, particularly the 
academic portion. And materials engineering 
per se has been left primarily to industry, which 
has yet to fully benefit from the expanded base 
of materials science. 

There is a twist in store. The science of 
materials has reached a level at which it now 
can radically change engineering practice. The 
possibilities are akin to what has come from 
the relationship between the life sciences and 
medicine. Until the 19th century, there was lit¬ 
tle or no science to guide medical technology 
and practice. Since then, however, the ever 
growing corpus of biomedical knowledge has been leading 
to an ever more amazing stream of health care innovations. 
Here, a genuine desire to meet societal needs has produced 
a healthy mix of reductionist and systems viewpoints, 
yielding a culture that naturally integrates scientific under¬ 
standing into practical use. The materials research commu¬ 
nity is poised to emulate this model. 


:Max PlahicicTor-.; 
miiUites the idea ; 
ptquante, therely 
s&jhg the stage ; 
fb^the develpp- 
ment of quantum 
mechanics. 
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Alfred Wilm dis¬ 
covers age hard¬ 
ening in alu¬ 
minum alloy, 
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used for making 
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other aircraft.. 
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Materials by Design: Efficient Innovation 

There is a general engineering design movement under way. 
It draws on the vast information pool generated by reduc¬ 
tionist analysis, but adds the component of design, for which 
the systems approach is crucial. 

.Central to the materials design approach is a powerful log¬ 
ical structure connecting the “four elements” of materials sci¬ 
ence: processing, structure, properties, and technological per¬ 
formance (see figure below). By connecting adjacent pairs of 
these elements, a three-link chain rq)resenting a versatile ma¬ 
terials paradigm emerges. The deductive, cause-effect logic of 
reductionist science flows from processing to performance. 
AH along the way, science reveals the relevant structures and 
phenomena, often in astounding clarity and detail. The induc¬ 
tive logic of systems engineering flows the other way, from 
performance to processing, thereby enabling designers to ar¬ 
rive at specific procedures likely to yield materials with the 
desired sets of properties and performance. 

Not often successful in terms of producing useful new ma¬ 
terials, early efforts at materials design nonetheless were 
harbingers of what materials development is to become. Gen¬ 
erally, what has made or broken past efforts was whether 
they included the eleihent of design. One notable suc¬ 
cess is the work of H.K.D.H. 

Bhadeshia and co-workers at 
Cambridge University, whose 
ambitious assault on the com¬ 
plex problem of weld metal 
design has spawned produc¬ 
tive efforts at several national 
laboratories in the United 
States and Japan. 

The initiative I know best 
is the Steel Research Group 
(SRG) at Northwestern Uni¬ 
versity, which my colleagues 
and I have continuously devel¬ 
oped since 1985 (5, 6). This university-industry-government 
program was organized within the context of systems engi¬ 
neering to explore general methods, tools, and databases for 
the design of materials, using high-performance steels as a 
test case (see sidebar on p. 995). There is a faith underlying 
this framework: The scientific knowledge base is now robust 
enough to supplant the traditional, empirically driven devel¬ 
opment of materials with a more efficient theory-driven and 
computationally based approach (see sidebar on p. 997). 

In the SRG, we begin by combining the perspectives of 
materials users, suj^liers, modelers, and designers into a set of 
specific materials propety objectives. Those specifications, in 
turn, help us define how to use, adapt, or expand science-based 
models and databases of material behavior. We then use these 
models and databases to zero in on compositions and process¬ 
ing protocols that can transform those compositions into alleys 
fitting our specifications. In the past decade, we have used this 
framework to develop new alloys with unprecedented sets of 
properties. Some are now under evaluation by industrial and 
government partners for use in airplanes, power generators, 
aircraft carriers, and other applications. 

Our own projects for designing steel are just a beginning. 
A recent NSF-sponsored workshop on Materials Design 
Science and Engineering (7) has called for broadening this 
approach to the design of all classes of materials. And in 
recognizing materials as one of five critical technologies for 
US. competitiveness, the President’s Office of Science and 
Technology Policy (8) has identified computational materi¬ 


als design as a principal opportunity. 

Opportunities abound for application of the new sys- 
tems-driven computational design approach. Successful ex¬ 
amples from our efforts include a stainless steel bearing for 
a space shuttle application, high-strength, high-toughness 
steels for aircraft landing gear and armor applications, and 
a new class of ultrahard steels for advanced gear and bear¬ 
ing applications. The more general validity of the design 
approach has been explored with computational design pro¬ 
jects focusing on nonferrous alloys, hydrate ceramics, case 
hardening polymers, and nanostructured thin film materials 
for microelectromechanical devices and hard coatings. 

The new design capabilities will also help realize the 
dream of biomimetic materials, which emulate the complex 
adaptive microstmetures of the living world that are beyond 
the reach of traditional empirical development. One success¬ 
ful demonstration is self-healing metalhc composites. These 
mcorporate “shape memory” alleys that exploit martensitic 
structural transformations to change their shape controUably. 
Integrating such components into thin film electronic devices 
to create efficient microactuators is leading to smart materials 
systems that unite the worlds of structural and elec¬ 
tronic materials. Here the growing philoso¬ 
phy of predictive materials design is now 
combining the electrical engineer’s realm of 
perfection-driven, artificially structured mi¬ 
crocircuits with the materials traditions of 
self-assembly (9) and defect tolerance (10). 

Following the same philosophy can re¬ 
duce the cost of discovery. In contrast to 
computational materials design, the preva¬ 
lent industrial methods of materials devel- 
opmoit are based on an mtrmsicaUy slow and expensive pro¬ 
cess of trial-and-error empiricism. Theoretical input has been 
qualitative at best. It typically takes tens of milhons of dollars 
over 2 decades to fully develop and quahfy a new material in 
a critical application (11). This sluggishness stands out in an 
era when engineers are expected to deliver new generations of 
products, such as automobiles, on an 18-month cycle. What’s 
more, under these competitive pressures many industrial 
materials-development groups have been severely reduced or 
closed over the past decade. A recent National Research 
Council study (11) has concluded that the greatest challenge 
to the materials field today is the short time and cost con¬ 
straints of the full materials development cycle. 

Another reason these industrial materials development ef¬ 
forts have been downsized is anchored in the old R&D mod¬ 
el in which new materials are discovered rather than de¬ 
signed. A typical estimate is that $1 of discovery costs $10 of 
development. If eliminating that first dollar was the only ad¬ 
vantage of computational materials design, its impact on the 
total discovery and development cost would be small. The 
real advantage of materials design is that good design in the 
first place requires much less development later. And rather 
than using materials design approaches simply to provide an 
initial prototype for subsequent empirical development, 
there’s yet more to be gained by integrating predictive model¬ 
ing throughout the full design and development process. 

Consider two costly and time-consuming phases of the 
standard materials development cycle: process optimization 
and qualification testing. A major concern in process optimiza¬ 
tion is scale-up. Because processing phenomena such as heat 
transfer dq)end on product size, a prototype material investigat¬ 
ed on a small scale is not likely to behave the same way when 
processed on a large scale. “Solidification design” can preempt 
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Model relay. Computational de¬ 
sign of materials requires a hier¬ 
archy of models. 


Virtualizing Materials to 
Create Real Ones 

Design of the hierarchical structure in materials 
requires a hierarchy of models based on materials 
science, applied mechanics, and even quantum 
physics (see figure at right). Models to help de¬ 
sign features on the coarsest structural level 
of solidification, such as the chemical iiicromechanics 
banding visible in the patterned sword 
blades mentioned earlier, employ 
powerful thermodynamic codes 
such as THERMOCALC (^3), 

These models enable de¬ 
signers to simulate the 
10-pm scale of structure. 

In metal alloys, this is the 
level at which the chemi¬ 
cal partitioning between 
liquid and solid phases 
evolves during solidifica¬ 
tion processing. Applica¬ 
tion of these models 
aids decision-making 
about thermal process¬ 
ing detaib, a practice we like to call "solidification design." 

At the 1.0-jim scale of structure, "transformation de¬ 
sign" is the goal. This concerns the evolution of structural 
changes during quenching, whereby crystal grains present 
at high temperatures transform and subdivide as the hot 
alloy cools into a hierarchy of lower temperature crys¬ 
talline units. The design goal here is to specify and control 
processing temperatures so that desirable microstructures 
will form, while hindering the formation of competing mi¬ 
crostructures that are less beneficial. 

The 0.1-^im scale represents the micromechanics de¬ 
sign level. An example of the phenomena relevant at this 
scale is "grain refining,” In which the large grains formed at 
high temperatures that can embrittle alloys are made 
smaller by more precise thermal or compositional control. 
There's a trade-off here, because more smaller particles 
can also catalyze ductile fracture (breaking, that is), as 
there are more interfaces that can separate from one an¬ 
other. Micromechanics models typically are based on 
continuum descriptions of mechanical phenomena 
that can simulate the evolution of microstructures 
during material deformation and fracture. 

In recent years, an even finer structural level—^the 
nanoscopic level—has become better understood and 
more controllable. The control of 1-nm-scale particle 
dispersions in alloys created through solid state pre¬ 
cipitation during "tempering” at intermediate temper¬ 
atures, for example, provides efficient obstacles for re¬ 
sisting plastic deformation. Said differently, this 
nanometer-scale structuring strengthens the metal. 

The development of design models for these 
diminutive scales builds on a half-century evolution 
of theory for both precipitation and strengthening in 
metals. A major materials science breakthrough of 
the 1950s was the identification of dislocations as 
the key defects that enable the sliding of crystal planes, 
which shows up as plastic deformation of many materials. 
While most structure-property relations, such as the Hall- 
Petch relation for grain refinement strengthening, are 


SolidificaHon 

Design 



1.0 


0.1 nm 


based on empirical 
correlations, a major 
triumph of the 1960s 
was the Orowan par- 
10 um tide strengthening 
equation. Derived di¬ 
rectly from disloca¬ 
tion theory, it relates 
strength to the inverse spac¬ 
ing of dislocation obstacles 
(such as nanoscale precipitates). 
Such developments in theory, which 
mathematically codify materials behav¬ 
ior, are cornerstones for the claim that 
materials can be designed largely in silico. 
The new accuracy of theoretical predic¬ 
tions of structure at the nanometer scale is 
only possible because of recent advances in high- 
resolution Instrumentation allowing precise calibration 
and validation of theory In this regime. This includes such 
techniques as x-ray and neutron diffraction, various electron 
microscopies, and atom-probe microanalysis. The latter is 
represented by the three-dimensional atomic reconstruction 
of a 3-nm strengthening carbide particle in an ultrahigh- 
strength steel shown below. Such new capabilities in struc¬ 
tural and chemical analysis down to the atomic scale open a 
new era of quantitative materials nanotechnology. 

The electronic level is the finest level relevant to real ma¬ 
terials. This is the realm of quantum design. As acknowl¬ 
edged by the 1998 Nobel Prize in chemistry shared by John 
Pople and Walter Kohn, the development of computational 
quantum mechanics and its extension via density functional 
theory constitute a profound advance that already has had 
significant industrial impact. A collaboration of materials sci¬ 
ence, applied mechanics, and quantum physics has enabled 
some of us in the Steel Research Croup (see main text) to 
apply computational quantum mechanics to engineer steel 
at the subatomic level. 

Our approach was to recast models of the impurity- 
induced embrittlement of grain boundaries into thermody¬ 
namic terms (14) and to rely on precise models of the 
atomic structures at grain boundaries. As represented by the 
computed valence charge density contours for a phosphorus 

atom at the core of an iron 
grain boundary, total energy 
calculations are sufficiently 
precise and accurate to ex¬ 
plain the known effects of 
Interstitial components such 
as boron, carbon, phospho¬ 
rus, sulfur, and hydrogen on 
the cohesion of iron grain 
boundaries. What's more, 
the calculations lead to new 
mechanistic insights at the 
electronic bonding level. Ex¬ 
tension to elements that 
occupy substitutional Fe 
sites in the boundary has 
enabled us to predict new 
alloying elements for enhancing boundary cohesion in 
steels. This is creating a new generation of "quantum steels” 
that incorporate these properties derived directly from elec¬ 
tronic-level predictions. -G.B.O. 



Atomscapes. A view of a Cr-Mo 
carbide particle that strengthens 
steel alloy. 
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this common showstopper. By using 
models that simulate how materials re¬ 
spond to processing conditions for de¬ 
signing a material at the ultimate process 
scale of interest, we can reduce the num¬ 
ber of expensive large-scale experiments 
while shortening the development cycle. 

Improving qualification testing is more challenging. To 
design with confidence, designers must know how properties 
vary within the materials they intend to use. Current practice 
requires much testing to define these variations statistically. 
But most structure-property theories provide mean values of 
the variables of interest; they do not tell a de¬ 
signer about location-to-location variations 
within materials. 

Help is coming, however. Researchers are 
developing a probabilistic materials science 
whereby structural distributions are mapped in¬ 
to property distributions. These are the kinds of 
data designers can use to make better predic¬ 
tions about how different materials will affect 
their products. Models of this kind already exist 
for fracture properties and heterogeneous phase transforma¬ 
tions, and more efforts ^e being planned for this vitally im¬ 
portant area. 

Structure of Education 

Education offers the most leverage for moving materials 
R&D into the systems-based paradigm. In its post-Cartesian 
form, modem metallurgy began with an emphasis on the di¬ 
rect correlation of processing and properties. The advent of 
physical metallurgy a century ago opened the “black box” of 
stmcture and brought a revolution in the fundamental under¬ 
standing of the mechanistic link between processing and 
properties. This understanding crea.ted the foundation for the 
more recent generalization to materials science, making pos¬ 
sible the general materials design and systems engineering 
methodologies described in this essay. Given the demonstrat¬ 
ed potential of materials design, institutionalizing it via edu¬ 
cation would result in many benefits to society. 

Due to the historic dominance of reductionist philosophy, 
however, we now have in place an analysis-oriented edu¬ 
cation system. In the name of objectivity, we tram 
students to shut down the more synthetic and emo¬ 
tional tools of thought, which are precisely the ones 
. best suited for doing good systems engineering. The 
future of engineering education, therefore, ought to tar¬ 
get these subjective reasoning powers. 

An activity of our new materials engineering curriculum 
at Northwestern, the Dragonslayer Project, is designed to do 
this. The project explores the design of an “aesthetic” mate¬ 
rial through collaboration of freshman and upperclass design 
teams. The widespread presence of Western dragon combat 
in literature is used to integrate frontier steel technology 
with the history and legend of swordmaking to design an ul- 
trahigh-performance sword that our market analysis shows 
would be of maximum value to sword collectors. 

Much of the mystique of the legendary S^urai sword, 
which even 5 centuries ago achieved performance levels 
equivalent to those of a modem carburized blade, stemmed 
from its ability to cut through other swords of the day. Die 
equivalent performance advance has been adopted as a 
benchmark for property requirements. 

Fanciful as the project might seem, it’s no different from the 
kmd of materials design that goes on at the DuPonts and Corn- 
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ings of the world. Searches of literature, 
albeit medieval hterature, were used to de¬ 
termine technical specifications, including 
the need for flame resistance when fight¬ 
ing fire-breathing dragons. Historians of 
the Dark Ages helped to identify and se¬ 
lect a historically accurate sword style— 
the patterned double-edge broadsword. And market analysis 
helped determine what the target buyers, sword collectors, 
would be willing to buy and how much they might pay. 

Finite-element mechanics simulations of cutting through 
a modem carburized blade led to a conceptual design for rad¬ 
ical surface-hardening technologies to achieve this ambitious 
sword-cutting objective. We also were guided by the unlikeli¬ 
hood that anyone would face supernatural evil armed only 
with technology. This led to an additional design parameter, 
which m^ket surveys indicated would constitute an impor¬ 
tant attribute in a collectible dragonslaying sword: The sword 
should be made from material of heavenly origin, namely, 
meteoritic iron. Diis is a technical feature with historical and 
legendary precedent (including the Excalibur legend). 

Our systems approach to materials design has estab¬ 
lished the feasibility of using tabletop-stale aqueous elec¬ 
trolytic refinement of available meteoritic iron to achieve 
the necessary purity for producing the unprecedented al¬ 
loy steel performance needed for a sword fit to defeat the 
most evil of adversaries. The freshman design team’s pro¬ 
posed market plan includes auctioning a single “techno- 
mystical” sword for publicity, followed by commercial 
marketing of a range of high-performance steel products, 
including the “Dragonslayer golf club.” 

In this millennium, a new architecture of synthetic 
thought will continue its symbiosis with modem computa¬ 
tional capabilities. And an age of empirical exploration will 
continue to be superceded by an Age of Design. This will 
open up powerful pathways by which human creativity, fused 
with scientific knowledge, will bring new levels of control 
over the material world widely applicable both to society’s 
problems and its ambitions. Die manifestation of the design 
paradigm, whether in an undergraduate or corporate setting, 
corresponds to a form of transmutation beyond the al¬ 
chemist’s dreams: the creation of materials from thought. 
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Matrix Alloys for Kinetic 
Energy Antitank Applications 
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Abstract- 


C 0 ^ f« iT: ' 

In recent years the US military has come under increasing eriti&isTTrregarding its environmental record. 
One area of contention has been the use of depleted uranium in its weapons and armors. As an antitank 
penetrator uranium-titanium (U- .075 wt. pet. Ti) has shown to be more effective thaycurrent tungsten (W) 
substitutes due to the density of the alloy and the effects of shear localization at the ^i^^^^^netration. It 
is the shear localization, which reduces the amount mushrooming in the head of the\^ap&Ti^ increasing the 
penetration depth. 


Despite several in-depth research studies no suitable material replacements have been found with this 
performance characteristic. However with current techniques, it is possible to design materials that not only 
undergo shear localization but that do so at the same point as the U-Ti alloy. 

Experiments were performed on Hf-Ta alloys to observe a martensitic to austenitic phase change. The 
atomic weight percent of the alloy was roughly calculated to get a To temperature for the martensitic to 
austenitic, a (HCP) to p (BCC), phase change, at 700°C using a straight-line method. Results from the 
experiments provide evidence that a martensitic phase structure is achievable using these calculations. The 
presence of this phase indicates that the subsequent transformation to austenite should be possible and that 
the system should provide the necessary shear localization properties required in a U-Ti replacement alloy. 

Further refinement of the Tq temperature calculation was performed with the use of ThermoCalc, a 
computer modeling program, providing a more accurate representation of the To temperature for the a 
(HCP) to P (BCC) phase change. Future research and experimental results will allow this model to become 
more accurate, to the point where a substitute alloy can be designed for the U-Ti system without any 
degradation in performance. 
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Figure 1: A cut a way view of a typical penetrator. Figure 2: The result of impact on W and U alloys (ARL) 

Introduction- 

While the end of the Cold War has reduced the threat of war, the battle on the materials 
science frontier rages on. The change in ro^ of the military as it becomes a policeman 
and peacekeeper, means that it must also change its approach to fighting. No longer is it 
acceptable to win at all costs. Increasing concern for the lives of its soldiers, civilians and 
the environment dictates the way in which the US military performs its current roles. 

Due to environmental concerns, the military has been trying to develop a tungsten (W) 
replacement for the depleted uranium-titanium (U-0.75 wt pet Ti) alloy that it currently 
uses in its kinetic anti-tank penetrators (Figure 1). Tungsten alloys, however, have shown 
to be considerably less effective in these applications, due to the plastic deformation that 
occurs when tungsten penetrators hit a target. The ability to design and model material 
systems, means that it is possible to replicate the properties of the U-Ti alloy in a 
tungsten system. 

The overall goal of the project is to be able to give a tungsten alloy the shearing 
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properties that will allow it to have similar or better performance characteristics than the 
existing U-Ti system. 

The approach is to design an alloy that exhibits shear localization (Figure 2) and that 
when used as a matrix alloy with tungsten will provide a system that exhibits overall 
shear localization. 

The matrix alloy is expected to compose around 30% of the tungsten penetrator. 

With this in mind the system that was chosen as the matrix alloy for these experiments 
was hafnium-tantalum (Hf-Ta) (Figure 3). The most important aspect of this system is the 
diffusionless phase change (a-P), which is expected to provide the shearing mechanism. 
Additional consideration was also given to the fact that these alloys share similar 
densities and processing characteristics with tungsten. 

The structure of the project is well outlined by Northwestern’s Materials Science 
department’s approach of looking at the relationship between Performance-Properties- 
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Figure 4: Project goals Figure 5; The stress vs. strain curve for U-Ti ^ 

Structure and Processing (Figure 4). The overall goal of the project was determined by 
the performance requirements set by the current U-Ti alloy. 

The major property of the U-Ti alloy system that is directly responsible for the current 
alloy’s overwhelming performance advantage is shear localization (Figure 5). This figure 
shows the stress verse strain curve for the actual U-Ti alloy, a typical W alloy, and a 
computer model of expected U-Ti curve. The dramatic drop in the actual stress as the 
stain increases means that excess metal shears off instead of plastic deforming and 
inhibiting further penetration.'^ 

The structure, which facilitates the shearing process, is a diffusionless phase 
transformation, in the case U-Ti p-y. In the Hf-Ta alloy system the phase change from a- 
P (BCC to HCP) is expected to result in a similar Effect. This Effect has been observed in 
steel systems where there is a martensitic to austenitic phase change. By manipulating the 
composition of the alloy the temperature at which phase change occurs can be selected. 
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The details of final subdivision, processing, were determined by^Hf-Ta system. The use 
of these refractory metals limits the options available to make a homogeneous alloy due 
to their high melting temperatures. This provides two possible options for the 
manufacturing of samples, powder sintering and vacuum arc melting. Powder sintering 
was used in these experiments. 

Literature Survey- 

Literature showed that quite a bit of previous work had been done both to analyze the U- 
Ti system and to develop tungsten substitutes. However at this time not much progress 
has been made in creating a tungsten alloy that meets the material performance 
requirements. The most important of these requirements being the ability of the 
penetrator alloy to undergo re-sharpening via shear localization. 

In order to simulate a ballistic/ impact the Hopkinson bar test appears to be favored due 
to the high compression rates that can be simulated. Testing on current U-Ti and W alloys 
has shown that this test produces accurate results relative to actual field performance. ^ 

No current tungsten heavy alloy, WHA, has been found that shows the shear localization 
effect. Mechanical testing of these alloy results in mushrooming at the head of the sample 
(Figure 2). U-Ti samples show a sudden failure and load drop in the stress strain curve, a 

V 

sign of the shear localization. This is due to the P to y phase transformation in the U-T/ 
system which occurs at a range in strain of .8 to 1.0 depending on the strain rate. Through 
mechanical testing it has been found that the strength of the U-Ti alloy drops significantly 




6 


when deformation proceeds in the high temperature y phase.® Current W alloys exhibit no 
such sinarl^ behavior. 

Research on steel has shown a correlation between a diffusionless martensitic phase 
(Image 1) change and early failure due to shear localization. The shear localization 
process relies on shear band formation and subsequent failure through shear band 
interaction^ductile rupture of the material. This causes the penetrator to efficiently eject 
excess material at the head of the weapon and penetrate further into armor material. 

Thermodynamically To temperature is the temperature at which the molar Gibbs energy 
of BCC and the HCP phases are equal. 

Gm (Ot)lxTa — Gjn (P)lxTa 

Above this temperature diffiisionless a to P transformation can take place. 

However partitioning by diffusion can still occur.’ Modeling the alloy from this 

/ 

thermodynamic information should provide the necessary phase change. 

Experimental Procedure- 

The samples were designed to achieve the martensitic to austenic phase change that is 
necessary for the shear localization. A rough approximation was made on the phase 
diagram using a straight-line method intersecting the point halfway between the Hfa and 
eutectic. This resulted in the following equation denoting the approximate To temperature 
for the phase change: 

y = -2803x+1743 
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Where x is weight percent tantalum and y is the To temperature in degrees Celsius. In 

■^siy' VxAf ■'i) tf ij- ' i £? ‘'^y 

order to achieve a To of 700°C the equation required 37.21 wt. pet. Ta. 

From this calculation Ta (13.9g) and Hf (29g) metal powders, -325 mesh (44|j,), were 
mixed using a ball milling process (Process flow chart 1). This was done in a N 2 
atmosphere with steel ball bearings for half an hour. The powder was then pressed using 
isothermal pressing at 50,000 lbs resulting in a sample that was approximately 20cm x 
3mm diameter. This sample was then divided into four smaller samples (12.2g, 9.9g, 
8.7g,6.4g) 


Two of these samples (12.2g, 9.9g) were then placed in a tungsten furnace and melted in 
an argon fpraaCe above 2600°C. This temperature was held for about five minutes 
(Process flow chart 2). These conditions were caused by a faulty reading from the 
thermocouple and caused the Hf-Ta system to form a liquid, essentially casting it in the 


graphite crucible. 

The formation of the alloy by this method caused some of the Ta components to separat^ 
creating an inhomogeneous mixture. This is further covered in the results and discussion 
section. Five samples, 10mm by 2mm in diameter, were cut from the cast piece using 
electrical discharge machining (EDM). These will be referred from here on out as the 


tungsten furnace samples. 



Image 2: Micrograph of sintered sample Image 3: SEM photo of pre-etched tungsten furnace sample 

surface 

Another attempt to sinter a sample was performed at 1850°C for Ihr. in an Ar atmosphere 
(Process flow chart 3). The physical constraints of the furnace required that the sample be 
smaller than the original size. This requirement resulted in a sample that was 1.5 cm long 
and 3mm in diameter. 

Since this sample was sintered for a relatively short time, Ihr. verse our target of 24hrs. 
the result of this was that the sample was not fully dense as the micrographs show 
(Image 2). 

Surface analysis of the sintered sample also showed several large cracks running through 
the sample. These likely occurred due to a temperature gradient inside the furnace. The 
porosity of the sample also made it very brittle and along with the surface cracks made 
processing for further testing difficult. This sample will be referred to as the sintered 


sample. 



Image 4: A post-etched SEM photo of tungsten Image 5; The microstructure after solution treatment 
furnace sample 

Energy Dispersive Spectroscopy was performed on a tungsten furnace sample using the 
3500 scanning electron microscope (SEM) at 20KeV both prior (Image 3) to and after 
(Image 4) etching. This was done to ensure that no contamination had occurred with the 
heat shields due to the high furnace temperatures. Counts were collected and analyzed by 
the Imix program (Process flow chart 4). Hf, Ta, Mo, C elemental K lines were inserted 
over the spectrum (Figure 5). 

The tungsten furnace samples were solution treated after some initial dilatometer testing. 
It was hoped that this would improve the microstructure. This was done in an argon 
atmosphere at 1200°C for Ihr. The samples were then quenched in oil (Image 5). 

One sample of each t 5 q)e (cast, cast-solutionized, sintered) was mounted in a cold mount 
and polished using 240, 320,400,600 grit sand paper followed by 30, 6,lp,m water based 
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Temp ®C 


Chart 1: Typical hysteresis from a phase change and a result from the tungsten furnace samples 


diamond polisher and .03|j,m alumina solution. Between each step the samples were 
cleaned with soap in an ultrasound and washed with methanol between polishing steps. 

To enhance visual analysis Hf etching of the grain boundaries was performed using a HF 
etchant (10ml HF, 45ml HNO 3 ,45ml H 2 O). The sample was swabbed for fifteen seconds 
to enhance the grain boundary definition. This assisted greatly in the analysis of the 
microstructure which became much more apparent after etching. Digital pictures were 
taken under 200,400, and lOOOx magnification. 


Three kinds of samples measurements were taken on the dilatometer, tungsten furnace 
pre-solution treated samples, tungsten furnace post solution treated samples, and a 
sintered sample. The dilatometer measures the relative change in length vs. the change in 
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temperature. It is from these experiments that a To temperature can be calculated from the 
hysteresis (Chart 1). By finding the austenitic start (As) temperature and martensitic start 
(Ms) temperature, denoted by a 1% change from linearity, the average of these two 
temperatures can be used to find an approximate To temperature. Further statistical 
manipulation can be performed to increase the calculations accuracy. 

A number of experiments were run in the dilatometer from 45°C to 1100°C on tungsten 
furnace samples, 10mm x 2mm, at the heating rate of 100°C/min. The samples were then 
held at 1100°C in an inert atmosphere for either 2 or 10 min., and rapidly or slowly 
cooled over a period of seconds (Table 1). An attempt to do a measurement of the 
sintered sample^ made. This sample was prepared by using t^ low speed saw to cut a 
wde off of the overall sample and then using the saw to ensure that the ends were parallel. 
However an accurate reading was not possible due to trouble attaching the thermocouple 
to the surface of the sample. 

Vickers hardness testing was performed on all three types of samples in the cold mount. 
The Vickers testing was performed under a 200g load for five seconds (Table 2). 

Results and Discussion- 

Both tungsten samples, pre and post solution treated, were unable to produce a hysteresis 
in the dilatometer. Despite the fact that the samples were solution treated to create a more 
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suspected Ta inclusions 


^ •' f j i i 

homogeneous mfetture^none of the trials showed enough deviation to determine an As and 
Ms temperature. Ten different trials were mn manipulating the heating rate, holding time, 
and cooling rate (Table 1). 

These results were due the inhomogeneous nature of the sample caused by phase 
separation due to the liquid to solid phase transformation, which produced Ta rich 
dendrites (Image 6). Also due to the short amount of time that the solution was in the 
liquid phase it is likely that not all the Ta powder metal melted into the liquid (Image 7). 
This was confirmed in the post solution treated sample which showed areas of 100% Ta 
during EDS analysis. Pulling Ta out of the system would have created a composition that 
had less wt. pet. Ta than indicated by the initial calculations making a diffusionless 
transformation more unlikely at 1100 °C 
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Atomic Fraction Ta 
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Figure 6: ThermoCalc Tq calculation for Hf-Ta ® Image 8: Sintered sample showing 

martensitic like structures 


Microstructural analysis of the post-solution treated samples shows some favorable 
results as far as the presence of a martensitic phase structure. When these photos are 
compared to that of martensitic regions in low carbon steel, the microstructure looks very 
similar (Image 7). This phase was not great enough to dominate the overall sample 
properties though since no phase change was detected in the dilatometer. 




The^p^dec^sintere^ample was very brittle and extremely porous due to the short 
sintering time that was dictated by the high temperature Ar furnace. Due to problems 
welding the thermocouple on the sintered samples no reliable dilatometer results were 
attainable. However microstructural analysis of the post-etched sintered sample also 
shows features that look martensitic (Image 8). 


While hardness data was difficult to ascertain and varied dramatically over the sintered 

^ ., r 

- / . 

specimen^db the porosity of the sample, Two microhardess tests showed a Vickers 
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Chart 2; A comparison of the straight-line approximation for To with the ThermoCalc modeling 


hardness of 753 and 841 significantly higher than the fully dense samples. This dramatic 

(7 ^ 'Y 

rise in hardness is alsojikel5rdue to a martensitic phase. 


The ThermoCalc To calculation from the Hf-Ta phase diagram results in the following 
equation (Figure 6). 

y = 2763.1x2 - 4897.4x + 1977.6 

Where x is the atomic fraction of tantalum and y is the To temperature in degrees Celsius. Y 
As expected this differs from the rough estimate provided by the straight-line method 
(Chart 2). Further refinement of this equation will be possible with experimental data. 
However this calculation provides considerable improvement as a starting point for future 
experiments. 
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Conclusions- 

The results of the testing on the Hf-Ta powder sintered sample show^ that it is possible 
to create martensitic regions in the system. This means that although a specific To 
temperature for the phase change was not confirmed experimentally this was a successful 
first step. The modeling of the To temperature from the Hf-Ta phase diagram allows for a 
more accurate composition to be calculated. Future work with the Hf-Ta system should 
allow for tailoring of a tungsten alloy to the penetrator application. 
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SEM (EDS) Analysis 



Pre-solution treatment EDS readings 


Element 

Point 1 

Point 2 

Point 3 

Point 4 

Average wt% 

Var 

stdev 

Ta 

0.00 

0.42 

0.64 

0.62 

0.42 

hee 


Hf 

llllillllllEEE 


0.13 


0.32 

■eee 





0.00 

0.00 

0.00 

0.00 

0.00 

C 

0.43 

0.24 

0.23 

0.14 

0.26 

0.01 

0.11 


Post solution treatment and etching EDS analysis 


Element 

Point 1 

Point 2 

Point 3 

Average wt% 

Var 

stdev 

Ta 


100.00 

100.00 

100.00 

0.00 

o.oo 

Hf 


0.00 

o.od 

0.00 

0.00 

0.00 

Mo 

0.00 

0.00 

o.oq 

0.00 

0.00 

0.00 


Note: EDS spectroscopy is not reliable for lighter weight elements when the accelerating voltage is greater 
than lOKeV. As a result carbon readings are inaccurate. ^ 
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Dilatometer Testing 


Trial # 

Sample 

# 

Soultionized 

Heating Rate 
°/min 

Temp 

°C 

Hold Time 
min 

Cooling Rate 
sec 

1 

1 

No 

100 

1100 

2 

2 

2 

2 

No 

100 

1100 

2 

2 

3 

3 

No 

100 

1100 

2 

2 

4 

4 

No 

100 

1100 

2 

2 

5 

5 

Yes 

100 

1100 

2 

2 

6 

1 

Yes 

100 

1100 

2 

120 

7 

2 

Yes 

100 

1100 

10 ^ 

120 

8 

3 i 

Yes 

100 

■n™ 

10 

120 

9 

4 

Yes 

100 n 

■lEEl 

2 

960 

10 

5 

Yes 

100 

1100 

10 

960 


Table 1 


Vickers Hardness Results 



testi 

test2 

tests 

test4 

Average 

Var 

stdev 

Sintered 

313 

166 

753 

841 

518.25 

81370.7 

285 

Tungsten 

Furnace 








Pre-Solution 

239 

273 

367 

435 

328.5 

5978.75 

77.3 

Post-Solution 

369 

303 

298 

369 

334.75 

1176.19 

34.3 


Table 2 
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Processing Flow Chart 1 and 2 


-Sample Preparation -Tungsten Furnace Treatment 
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Processing Flow Chart 3 and 4 

-High Temperature Ar Furnace Treatment -Testing Routine 


























